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Traditionally models of accretion of gas on to T Tauri stars have assumed a dipole stellar magnetosphere, partly for sim- 
plicity, but also due to the lack of information about their true magnetic field topologies. Before and since the first magnetic 
maps of an accreting T Tauri star were published in 2007 a new generation of magnetospheric accretion models have been 
developed that incorporate multipole magnetic fields. Three-dimensional models of the large-scale stellar magnetosphere 
with an observed degree of complexity have been produced via numerical field extrapolation from observationally derived 
T Tauri magnetic maps. Likewise, analytic and magnetohydrodynamic models with multipolar stellar magnetic fields have 
been produced. In this conference review article we compare and contrast the numerical field extrapolation and analytic 
approaches, and argue that the large-scale magnetospheres of some (but not all) accreting T Tauri stars can be well de- 
scribed by tilted dipole plus tilted octupole field components. We further argue that the longitudinal field curve, whether 
derived from accretion related emission lines, or from photospheric absorption lines, provides poor constrains on the 
large-scale magnetic field topology and that detailed modeling of the rotationally modulated Stokes V signal is required 
to recover the true field complexity. We conclude by examining the advantages, disadvantages and limitations of both the 
field extrapolation and analytic approaches, and also those of magnetohydrodynamic models. 
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1 Introduction 

The original models of magnetospheric accretion on to a T 
On Tauri star assumed the stellar magnetosphere was a dipole 
^ (Konigl 1991; Collier Cameron & Campbell 1993; Hart- 
mann et al 1994; Shu et al 1994). This assumption was made 
. partly for simplicity and partly due to a lack of available in- 
■ formation about the true large-scale magnetic field topology 
• . of such stars. None-the-less dipole magnetospheric mod- 
^ ^ els have proved successful in reproducing many observa- 
tional signatures, for example, the profiles and absorption 
^ ' components of accretion related emission lines (e.g. Hart- 
mann et al 1994; Kurosawa et al 2006). Since the middle of 
the 1990s several two dimensional magnetohydrodynamic 
(MHD) models have been presented (e.g. Goodson et al 
1997; Miller & Stone 1997; Romanova et al 2002; Kiiker 
et al 2003; von Rekowski & Piskunov 2006; Bessolaz et al 
2008; Zanni & Ferreira 2009), as well as some three dimen- 
sional (3D) simulations with tilted dipole magnetospheres 
(Romanova et al 2003, 2004; Orlando et al 201 1) - see Gre- 
gory et al (2010) for a review. 

Dipolar accretion models are not without their critiques 
(e.g. Safier 1998), and dropping the dipole assumption can 
significantly affect the structure of accretion columns within 
the magnetosphere (Gregory et al 2006a), the locations of 
hot spots at the base of accretion funnels (Gregory et al 
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2005; Mohanty & Shu 2008), and the density of the gas in 
the flow (Gregory et al 2007). High order field components 
may even play a dominate role in the physics of the gas in- 
flow as the accretion columns approach the star (Adams & 
Gregory 2011). 

Through the analysis of Zeeman broadened lines in in- 
tensity spectra T Tauri stars have been found to host surface 
averaged magnetic fields of order a kilo-Gauss (see Johns- 
Krull 2007 and references there-in). Provided that the sur- 
face fields are associated with a stellar magnetosphere that is 
sufficiently globally ordered then such strong fields are able 
to disrupt the circumstellar disk at a distance of a few stellar 
radii, as required by the magnetospheric accretion scenario 
(Konigl 1991). There appears to be little variation in the 
average surface field between sources, and between accret- 
ing and non-accreting stars, although intriguingly the aver- 
age (unsigned) magnetic flux appears to be less, on average, 
for stars in older star forming regions (Yang & Johns-Krull 
2011). 

Spectropolarimetry provides a method of probing the 
magnetic field topology. By measuring the circular polariza- 
tion signal in magnetically sensitive photospheric absorp- 
tion lines, and in accretion related emission lines, informa- 
tion can be obtained about the surface fields of T Tauri stars. 
Zero net polarization signals have often been measured in 
photospheric absorption lines (e.g. Valenti & Johns-Krull 
2004). As such lines form essentially uniformly across the 
entire visible stellar surface, this suggests that T Tauri stars 
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host complex surface fields with roughly equal amounts of 
positive and negative field regions. However, often accre- 
tion related emission lines, such as Hel 5876A, show strong 
levels of circular polarization (e.g. Johns-Krull et al 1999), 
with the signal modulated on a timescale equal to the photo- 
metrically determined rotation period of the star (Valenti & 
Johns-Krull 2004). This suggests that for many T Tauri stars 
magnetospheric accretion occurs predominantly into strong 
single polarity spots on the surface the star. However, such 
studies which are based on measuring the phase modulated 
stellar disk averaged longitudinal field component, cannot 
constrain the field topology of the star beyond the interpre- 
tation that T Tauri stars have complex surface fields, and 
somewhat simpler large scale fields. 

In 2007 the first magnetic maps of the large-scale field 
of an accreting T Tauri star were published (Donati et al 
2007). At the time of writing magnetic maps have now been 
produced for 10 accreting T Tauri stars namely V2129 Oph, 
BP Tau, V2247 Oph, AA Tau, TW Hya, both stars of the 
close binary V4046 Sgr, CR Cha, CV Cha and MT Ori (Do- 
nati et al 2007, 2008b, 2010a,b, 2011a,b,c; Hussain et al 
2009; Skelly et al 2011). They are produced by monitoring 
the rotational modulation of the circular polarization sig- 
nal over a complete stellar rotation cycle, and in practise 
several, detected in both photospheric absorption and accre- 
tion related emission lines. As the polarization signal in in- 
dividual photospheric absorption lines is vanishingly small 
cross-correlation techniques (such as Least-Squares Decon- 
volution (LSD); Donati et al 1997) are employed in order 
to extract information from as many of the spectral lines as 
possible. The complete process, known as Zeeman-Doppler 
imaging, see Donati (2001) for a review and Donati et al 
(2010b) for its specific application to accreting T Tauri stars, 
simultaneously reconstructs maps of the radial, azimuthal 
and meridional field components at the stellar surface, as 
well as a brightness map (the surface distribution of cool 
spots) and a map of excess emission in accretion related 
lines (the surface distribution of accretion hot spots). Due 
to the flux cancellation effect that affects all polarization 
measurements, the magnetic maps only contain informa- 
tion about the large scale field topology. Once the magnetic 
maps have been derived the field components can be written 
in terms of spherical harmonics (e.g. Donati et al 2006b). 
From the maps the coefficients of the spherical harmonic 
decomposition can be calculated, and from these the polar 
strength of the each of the multipole moments along with 
their tilts relative to the stellar rotation axis, and the phase 
which they are tilted towards, can be determined. A descrip- 
tion of the observational process involved in constructing 
stellar magnetic maps via Zeeman-Doppler imaging, as well 
as the assumptions and limitations, can be found elsewhere 
in this volume in the paper of G. A. J. Hussain, to which 
readers are referred for further details (see also Donati & 
Landstreet 2009). 

Most of the published magnetic maps have been ob- 
tained as part of the Magnetic Protostars & Planets (MaPP) 



program. This ongoing large program consists of 690 hours 
of observing with the ESPaDOnS (Echelle SpectroPolari- 
metric Device for the Observation of Stars; Donati et al 
2006a) instrument at the Canada-France-Hawai'i Telescope 
(CFHT) between 2008 and 2012, 1 130 hours with the NAR- 
VAL spectropolarimeter (Auriere 2003) at the Telescope 
Bernard Lyot (TBL) in the Pyrenees, and additional ground 
and space based observations (e.g. XMM-Newton, Chandra 
and HST) for select sources. 

In this paper we discuss the results from the ongoing 
MaPP program and present magnetic maps for some accret- 
ing T Tauri stars in section [21 We demonstrate that the large 
scale field of some (but not all) T Tauri stars are well de- 
scribed by fields consisting of a tilted dipole plus a tilted 
octupole component. We review the current 3D field extrap- 
olation and analytic models of their magnetospheres in sec- 
tion [3l and argue that use of the stellar disk averaged lon- 
gitudinal field component alone provides poor constrains 
on the magnetic field topology. We conclude in section |4] 
by comparing both approaches to one another, and also to 
magnetohydrodynamic models, discussing the advantages, 
disadvantages and limitations of each. 

2 Surface magnetic maps 

T Tauri stars are found to host multipolar magnetic fields, 
but the field topology varies between stars and over time for 
the same star. All appear to have dipole components that 
are strong enough to disrupt their circumstellar disks. Al- 
though their large scale fields are non-dipolar, they are often 
simpler than the highly complex fields of rapidly rotating 
zero-age main sequence stars (e.g. Donati et al 2003). This 
result is not a phase coverage effect, as demonstrated by 
Hussain et al (2009). Some T Tauri stars stars host simple 
axisymmetric fields with strong kilo-Gauss dipole compo- 
nents (e.g AA Tau and BP Tau). Others host dominantly ax- 
isymmetric fields but field modes of higher order than the 
dipole are found to be the most dominant (typically, but 
not exclusively, the octupole component), and their dipole 
components may be strong or weak (e.g. TW Hya, V2129 
Oph and MT Ori). Others host complex magnetic fields that 
are highly non-axisymmetric with very weak dipole com- 
ponents «- 0.1 kG; e.g. V4046 Sgr AB, CR Cha and 
CV Cha). There is strong link with stellar internal struc- 
ture where fully convective stars close to the end of the 
fully convective phase host simple fields with strong dipole 
components. The dipole component decays rapidly with the 
development of a radiative core, although the large scale 
field appears to retain its axisymmetry until the core oc- 
cupies a substantial volume of the star. Empirically once 
the mass of the core exceeds about 40% of the stellar mass, 
the field becomes highly complex and the dipole compo- 
nent is weak (see Gregory et al 2012 for full details of how 
T Tauri magnetic topologies vary with stellar internal struc- 
ture). There is also tentative evidence for a bi-stable dynamo 
process operating amongst the lowest mass fully convective 
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Fig. 1 Magnetic maps showing only the radial component of the field of four accreting T Tauri stars derived from 
Zeeman-Doppler imaging. The large-scale fields of AA Tau in January 2009 (upper left; Donati et al 2010b), BP Tau in 
February 2006 (upper right; Donati et al 2008), V2129 Oph in July 2009 (lower left; Donati et al 2011a) and TW Hya in 
March 2010 (lower right; Donati et al 2011b) are multipolar, but well described as dipole-octupole composite fields with 
different polar strengths and moment tilts for each star (see Table [T]). The maps are shown in flattened polar projection with 
positive/negative field regions in red/blue and labeled in Gauss. The stellar rotation pole is in the middle of each map, the 
dotted lines are lines of constant latitude separated by 30° with the equator as the bold circle. The numbers and tick marks 
around the circumference denote phase and phases of observation respectively. Other T Tauri stars possess more complex 
large-scale magnetic fields, see Figure [2l 

Table 1 Polar field strengths of the dipole and octupole components, their tilts relative to the visible stellar rotation pole 
(/3), and the rotation phase that they are tilted towards (<l>) for the stars show in Fig.[T] The values listed for BP Tau were 
derived using an experimental version of the magnetic imaging code and have recently been revised from the values shown 
here (see footnote 1 for details). The fields of TW Hya and V2129 Oph were significantly different at other observing 
epochs. 
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Pdip 


^dip 
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AATau 


Jan09 


1.7 


170° 


0.65 


0.5 


10° 


0.50 


0.29 


BP Tau 


Feb06 


1.2 


20° 


0.65 


1.6 


10° 


0.15 


1.33 


V2129 Oph 


Jul09 


1.0 


15° 


0.40 


2.2 


20° 


0.50 


2.20 


TW Hya 


Mario 


0.7 


5° 


0.25 


2.8 


>175° 


0.85 


4.00 



T Tauri stars, where we expect some to host complex fields, 
and others simple fields with strong dipole components like 
the more massive fully convective pre-main sequence stars. 
This behavior mirrors the magnetic topology trends already 
found for main sequence M-dwarfs on either side of the 
fully convective limit (Donati et al 2008a; Morin et al 2008, 
2010, 2011). 

Figure [T] shows four examples of magnetic maps of ac- 
creting T Tauri stars. Only the radial field component is 



shown, and we have deliberately selected those stars which 
have large-scale field topologies that are somewhat simple 
and well described as dipole-octupole fields, as these will 
be useful for our analytic work in section 13.21 Other stars 
have significantly more complex fields, see Figured Maps 
of the other field components along with brightness and ex- 
cess accretion related emission maps can be found in the 
appropriate papers. 
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Fig. 2 As Figure[T]but showing magnetic maps of accreting T Tauri stars with large scale fields that are not well described 
as dipole-octupole composite fields; V2247 Oph in July 2008 (upper right; Donati et al 2010a), MT Ori in December 2008 
(upper left; Skelly et al 2011), and both stars of the close binary V4046 Sgr in September 2009 (Donati et al 2011c), the 
primary (lower left) and the secondary (lower right). The field of MT Ori is similar to that of TW Hya and V2129 Oph in 
that its field is dominantly axisymmetric but field modes of higher order than the dipole dominate (in this case there are 
roughly equal contributions from the octupole [£ = 3], the dotriacontapole [£ = 5] and the £ = 7 field modes). The fields of 
the other stars are complex and non-axisymmetric with weak dipole components, as are the fields of CR Cha and CV Cha 
(Hussain et al 2009; maps not shown here). 



The four stars shown in Figure [T] have dipole and oc- 
tupole components of differing polar strength, and tilts rel- 
ative to the stellar rotation axis, as summarized in TablellQ 
In this work we define the tilt (3 of the dipole or octupole 
moment to be that of the main positive pole relative to the 
visible rotation pole of the star. Therefore a star which has 
the main negative pole of the dipole coincident with the vis- 
ible rotation pole would have Pdip = 180°, and likewise 
Poet = 180° would correspond to the main negative pole of 
the octupole coinciding with the visible rotation pole. This 
is applicable to AA Tau and TW Hya which are close to a 
configuration with anti-parallel dipole and octupole compo- 
nent at the epochs listed in Table [T] An alternative notation 
would be to define the tilts of the components according 
to whichever main magnetic pole was in the visible hemi- 

^ The values listed in Table[T]for BP Tau were derived before the MaPP 
program began using an experimental version of the magnetic imaging 
code. A preliminary reanalysis of this older BP Tau data using a more 
mature version of the code (to presented alongside newly obtained data 
and new magnetic maps) suggests that B^ip = 1.3 kG, ^dip = 5°, 
^f^ip = 0.85, Boct = 1.8 kG, ^oct = 5° and ^oct = 0.15 are more 
appropriate for this star in February 2006. As the new data has not yet been 
published we have used the old values and magnetic map for BP Tau in this 
paper. 



sphere of the star, and then list the polar strength of that 
component as negative. For example, using this alternative 
approach for AA Tau then Bdip = — l.TkG, /3dip = 10° 
and ^dip = 0.15 (since tilting the negative pole of the dipole 
by 10° towards phase 0.15 is equivalent to tilting the posi- 
tive pole by 170° towards phase 0.15 + 0.5 = 0.65). Con- 
servative estimates of the uncertainties in the tilts (/3) and 
the phases of the tilts (^) listed in Table [T] are ^ 10° and 
0.1 (especially for moments with small tilts relative to the 
stellar rotation axis) respectively. 

All four stars in Figure [T] and Table [T] have been ob- 
served at two different epochs. V2129 Oph showed signif- 
icant changes in the strength of the dipole and octupole 
field components between observing epochs with the dipole 
component increasing by a factor of about three between 
2005 and 2011 (from 0.3 kG to 0.9 kG, see Donati et al 
201 la). Likewise for TW Hya where the dipole component 
appeared to flip from roughly anti-aligned with the octupole 
component at one epoch, to being titled by about 45° from 
an anti-parallel configuration at the other epoch; although 
as noted in Donati et al (201 lb) this may be due to the poor 
phase coverage at one epoch and requires additional obser- 
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vations to determine if such large scale topology changes 
are genuine and common for TW Hya. AA Tau and BP Tau 
showed little change in their field topology between observ- 
ing epochs, save an overall quarter phase shift in the entire 
map of BP Tau. This was most likely due to a small error 
in the assumed rotation period building up over the roughly 
40 stellar rotations between observing epochs. Changes in 
the field topology over time for both stars cannot be ruled 
out and further observations are required. We note that the 
magnetic maps of BP Tau were constructed and published 
using an experimental version of the magnetic imaging code 
for accreting T Tauri stars. The Zeeman-Doppler imaging 
process for such stars involves the use of polarization in- 
formation in accretion related emission lines, and therefore 
differs from the magnetic imaging process as it applies to 
other stars which lack such proxies. 

3 Multipolar magnetospheric models 

In anticipation of the newly available data stream from the 
ESPaDOnS spectropolarimeter and the observationally de- 
rived magnetic maps of accreting pre-main sequence stars, 
Gregory et al. (2005) produced the first simulations of the 
magnetospheres of T Tauri stars with complex magnetic 
fields. The magnetic maps of young rapidly rotating zero- 
age main sequence stars were used to generate 3D field 
extrapolations which where then surrounded by accretion 
disks. Consequently Jardine et al (2006) and Gregory et 
al (2006a,b) demonstrated that magnetic fields with an ob- 
served degree of complexity were required to explain many 
observational results, from the increase in X-ray luminosity 
with stellar mass, to the common detection of rotationally 
modulated X-ray emission (Flaccomio et al 2005). Subse- 
quently 3D MHD simulations of the star-disk interaction 
with non-dipolar fields were published by Long et al (2007, 
2008) using dipole-quadrupole fields, by Long et al (2009) 
using dipole-octupole fields, and by Romano va et al (201 1) 
and Long et al (201 1) using dipole-octupole fields with po- 
lar field strengths and tilts that match those derived from 
Zeeman-Doppler imaging studies of the accreting T Tauri 
stars V2129 Oph and BP Tau. In addition to the field ex- 
trapolation models, and the MHD simulations, multipolar 
magnetic fields can also be handled analytically (or at least 
semi-analytically), which allows for greater clarity of un- 
derstanding of the physical results from numerical models 
(Gregory et al 2010). Mohanty & Shu (2008) have also pre- 
sented a semi-analytic generalization of the Shu X-wind 
model which drops the dipole magnetosphere assumption. 

The spectropolarimetric results, and the derived mag- 
netic maps, provide clear motivation for constructing 3D 
models of the non-dipolar magnetospheres of pre-main se- 
quence stars. In this section we provide a detailed review of 
two of the three approaches to handling multipolar stellar 
magnetic fields, namely field extrapolations from magnetic 
maps and analytic models (we compare these to MHD mod- 
els in section 131). Although our attention is focussed on the 



magnetic fields of accreting pre-main sequence stars the re- 
sults that we derive and discuss herein are applicable gener- 
ally to the magnetosphere of any star or planet. 

3.1 Field extrapolations from magnetic maps 

Stellar magnetic maps derived from Zeeman-Doppler imag- 
ing can be used as inputs to models of the 3D large-scale 
magnetosphere. Numerical field extrapolations from mag- 
netic maps have been produced for stars of various spectral 
type and evolutionary stage, for example, zero-age main se- 
quence rapid rotators (e.g. Jardine et al 2002); accreting T 
Tauri stars (Jardine et al 2008; Gregory et al 2008); solar- 
type post T Tauri stars (Dunstone et al 2008; Marsden et al 
2011); a massive B star (Donati et al 2006b); and for the 
global solar magnetic field (e.g. Riley et al 2006; Ruan et 
al 2008). The potential field source surface method is used 
to produce the field extrapolations, a technique originally 
developed by the solar physics community (Altschuler & 
Newkirk 1969; Schatten et al 1969). A potential field repre- 
sents the lowest possible energy state of the stellar magnetic 
field. Twisting of the field line by stellar surface transport 
effects (surface differential rotation, meridional flow etc) 
leads to a departure from the potential state and allows mag- 
netic energy to be stored in the field which can eventually be 
released through magnetic reconnection events that trigger 
flares. 

The potential field source surface model has two bound- 
ary conditions. The first is that the radial field component at 
the stellar surface is the same as determined from the obser- 
vationally derived magnetic map of the star whose magne- 
tosphere is being modeled (the magnetic map is thus used 
as a direct input). The second boundary condition is that the 
field is purely radial (Bq = 5^ = 0) at a spherical equipo- 
tential surface at some height above the star. This surface, 
of radius Rs, is known as the source surface and provides 
a simple method of mimicking the stellar wind plasma dis- 
torting and pulling open the closed field beyond a certain 
height above the star|3 The source surface thus represents 
the maximum extent that a closed field line loop can have, 
and adjustment of this boundary affects the amount of open 
and closed flux through the stellar surface reconstructed in 
the field extrapolation model. Smaller (larger) values of i^^ 
yield more (less) open field relative to closed field (see fig- 
ure 1 of Gregory 2011). 

In the solar case the source surface radius is constrained 
by in-situ heliospheric satellite observations to be Rs ~ 
2.5 Rq. For stars such in-situ satellite observations are not 
available and other indirect methods must be used to con- 
strain the radius of the source surface. Although various 
ideas and suggestions have been put forward, Rs remains 

^ Potential field source surface models have been produced with non- 
spherical source surfaces (e.g. Schulz et al 1978). In the case of the Sun 
the source surface is a prolate spheroid at solar minimum with major axis 
along the solar rotation axis. At solar maximum the source surface is more 
complex but is approximately spherical once averaged over polar and az- 
imuthal angles (Riley et al 2006). 
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a poorly constrained parameter for stars other than the Sun. 
However, once a value for Rs has been chosen, a coronal 
X-ray emission model can be used to calculate the global 
X-ray emission measure for comparison with independently 
obtained X-ray observations (Jardine et al 2006; Hussain et 
al 2007). Adjustments in Rs affect the amount of closed 
field, and consequently the amount of X-ray bright regions 
on the star. However, once Rs exceeds a few stellar radii the 
change in the calculated X-ray emission measure (or equiv- 
alently the X-ray luminosity) is small for successively larger 
Rs values (Jardine et al 2008). For accreting T Tauri stars 
surface hot spots, which form at the shocks at the base of 
accretion columns, can be used as an additonal constraint 
on Rs. Maps of the excess emission in certain accretion re- 
lated emission lines are constructed as part of the Zeeman- 
Doppler imaging process (see Donati et al 2010b). Current 
observations suggest that for the majority of accreting T 
Tauri stars hot spots are at high latitude towards the visi- 
ble rotation pole (see also Petrov et al 2011). In order to 
explain how material can accrete from the disk, in the stel- 
lar midplane, into high latitude hot spots the source surface 
must be large enough to ensure that closed field can reach 
towards the polar regions. 

Potential fields are current free and therefore at every 
point exterior to the star V x B = 0. Assuming that the field 
external to the star can be written as the gradient of a mag- 
netostatic scalar potential ^, B = — V^, then this com- 
bined with Maxwell's equation that the field is solenoidal 
(V • B = 0) means that the magnetostatic potential at any 
point (r, 6>, (j)) exterior to the star can be determined by solv- 
ing Laplace's equation V^^ = 0. Thus the spherical field 
components (e.g. Br = —d^/dr etc) can be calculated at 
each point exterior to the star, and then a field line tracing 
algorithm employed to reconstruct the 3D global magnetic 
field of the star. Jardine et al (2002) provide further details 
of the potential field source surface method as applied to 
stellar magnetic maps, with Jardine et al (2006) and Gre- 
gory et al (2006a, 2008) providing details of its application 
specifically to T Tauri stars. Potential field models produce 
unique 3D field topologies (Aly 1987), once a source sur- 
face radius has been selected, however the fields are static 
in time. Such models cannot be used to model the time evo- 
lution of the magnetic field which will evolve due to motion 
of the foot points on the stellar surface (surface transport 
effects) and by the twisting of the large field due to the in- 
teraction with the disk. We return to these points in section 

m 

Figure [3] shows potential field extrapolations from the 
observationally derived magnetic map of the accreting T 
Tauri star V2247 Oph (Donati et al 2010a). This star has 
a complex magnetic field (see the map in Figure [J]) with 
a weak dipole component (< 0.1 kG); none-the-less the 
dipole component dominates on the large scale (but see sec- 
tion [3]2]2|). A source surface of i^^ ~ 4.3 i?*, where R^ = 
1.6 R©, comparable to the equatorial corotation radius for 
V2247 Oph (Rco = 4.25 i?* with = 0.36 M©, Prot = 

© 201 1 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 
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Fig. 3 Field extrapolations from the magnetic map of 
V2247 Oph, shown in Figure [2] and the bottom row, show- 
ing the smaller scale field that would contain high temper- 
ature plasma that constitutes the star's X-ray bright corona 
(top row), the larger scale field that interacts with the inner 
disk funneling gas on to the stellar surface (second row), 
and open field lines along which a stellar wind could be 
launched (third row). The star is viewed at inclination 45° 
and at phase 0.5. The top image has been expanded for clar- 
ity and is not on the same scale as the others. The white 
points in the plot of latitude versus longitude are the foot- 
points of the accretion funnels showing where the hot spots 
would be on the stellar surface, assuming an anti- symmetric 
field distribution with respect to the stellar midplane (see 
section [3i2JJ . www.an-journal.org 
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3.5 d) has been chosen for illustrative purposes. The radius 
is different from that in Donati et al (2010a) as we have 
updated the distance estimate to V2247 Oph from 140 pc 
to 120 pc (Loinard et al 2008), which yields a different lu- 
minosity and consequently a different radius. Assuming ac- 
cretion occurs from a range of radii within corotation the 
bottom row of Figure [3] shows the location of the accre- 
tion spots on the stellar surface. The main spot in the vis- 
ible hemisphere is arc- shaped and decreases in latitude to- 
wards the stellar equator from phase ~ 0.4 to ~ 0.8. This is 
in excellent agreement with the accretion spot location in- 
ferred from the excess emission map derived by Donati et al 
(2010a; see their figure 7). 

3.2 Analytic models 

Multipolar magnetic fields can also be handled analytically, 
including highly complex magnetic fields with many high 
order components. For the analysis in this work we limit 
ourselves to the zonal (also called axial) multipoles of arbi- 
trary I number (or linear combinations of I numbers). Thus 
we are neglecting the non-axisymmetric field modes (those 
with m 7^ 0). The same assumption is made in the gener- 
alized multipolar X-wind model of Mohanty & Shu (2008) 
and in the initial field structures used as inputs in current 3D 
MHD simulations (Romano va et al 2011; Long et al 2011; 
see section |4]). Observationally the magnetic fields of ac- 
creting T Tauri stars are dominantly axisymmetric with the 
exception of some low mass stars (typically below 0.5 M©) 
and older/more massive T Tauri stars with substantial ra- 
diative cores (those with Mcore/^^ ^ 0.4; Gregory et al 
2012). 

Gregory et al (2010) derived an expression for the spher- 
ical magnetic field components {Br^ Bq^ B^j)) exterior to a 
star for a multipole of order £ (where ^ = 1, 2, 3, 4, ... are 
the dipole, the quadrupole, the octupole, the hexadecapole 
etc). 



B = 



■^i,pole 
N 



£+2 



k=0 



(_l)fc(2£_2fc)! 
2'''k\{e - k)\{e - 2k)\ 



-2k 



[{2£ -2k^ l)r -{£- 2k){fi • r)" V] }, 



(1) 



where N = £/2 or (^-1) /2 whichever is an integer, Bi'^''^'' 
is the polar field strength of the ^-number multipole, f is a 
unit vector along the radial direction and /x is a unit vector 
in the direction of the multipole moment fj.. Magnetic fields 
consisting of multipoles with many high order ^-number 
components can be constructed by summing equation ([T]) 
over the £ numbers of interest. 

In the following subsections we derive expressions for 
the spherical field components that are written in terms of 
the tilt and phase of tilt of each multipole moment explic- 
itly. Using these expressions we construct simulated mag- 
netic maps for the four stars shown in Figure[Tl all of which 



have large scale fields that are well described by tilted dipole 
plus tilted octupole components. We compare the simulated 
maps to the maps derived from Zeeman-Doppler imaging 
which contain information about many higher order £ num- 
ber and m 7^ field modes. We further derive generalized 
expressions for the stellar disk averaged longitudinal field 
component and demonstrate that reliance on this sole diag- 
nostic provides poor constrains on the overall surface mag- 
netic field topology (e.g. Donati & Landstreet 2009). 

3.2.1 Tilted multipolar magnetic fields 

We consider an axial multipole (m = 0) of order £ tilted by 
an angle /3 relative to the stellar rotation axis and tilted to- 
wards a longitude i/j where is related to the rotation phase 
^ via 

= (1 - ^) X 360°. (2) 

Note that rotation phase $ and longitude i/j run backwards 
such that longitudes = 0°, 90°, 180°, 270° and 360° cor- 
respond to rotation phases ^ = 1, 0.75, 0.5, 0.25 and zero, 
and that our derivation assumes the xz-plane corresponds 
to a rotation phase of zero (see Appendix lAb. The spheri- 
cal field components can be derived from equation O, fol- 
lowing the details of Appendix [Al and by noting that as 
fi = ii~^{iir,l~ie,IJ^cj?j then/i • r = /ir/^, 

£+2 



Br = 5f '^^^"^ 



N 



{-lf{2£-2k)\ 



f^^y2^k\{£-k)\{£-2k)\ \ii 



(3) 



Be 



B^ 



i,pole / D \ ^+2 



^ [ {-lf^\2£-2k)\ 



k=0 



Bfh = 



2^k\{£ - k)\{£ - 2k - l)\ \ii 

1^^ I -tt^ 



-2k-l 




N 



2. 



e-2k-l 



^ .(5) 



where we again note that N = £/2 or {£ — l)/2 whichever 
is an integer, /i = and the components of the multipole 
moment are 

jiir = 1^ sin cos (j) cos ip sin /3 + /i sin sin cj) sin sin P 
+ /i cos cos P (6) 

/jj0 = fi cos 6 cos (j) cos i/j sin /3 -\- fi cos 6 sin (j) sin i/j sin /3 
— ji sin cos j3 (7) 

/i0 = — /i sine/) cos sin /3 + /i cos (/) sin sin /3 (8) 

where and (j) are standard spherical coordinates in the ref- 
erence frame of the star. The components of a multipolar 
magnetosphere can then be determined by summing the ap- 
propriate field components. For example, a dipole plus oc- 
tupole magnetosphere where the dipole (octupole) moment 
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of polar strength Bdip (Boot) is titled by /Sdip {Poet), relative 
to the stellar rotation axis, towards a rotation phase denoted 
by V^dip i^poct) has spherical components 

+ Br{BocUpocU^ocU^ = ^) (9) 

Be = Bo{Bdip, Pdip.i^dip,^ = 

+ Be{BocuPocu^octJ = 3) (10) 

B(j) = B(f){Bdip^ Pdip^i^dip^^ = ^) 

+ B^{Boct.PocU^octJ = 3) (11) 

where Bd^p = bI^^'^' and Boat = bI^^'^\ 

We note that the field components derived in this section 
for tilted axial multipoles (those with m = 0) also apply 
to sectoral multipoles (those with i = m) which are axial 
(zonal) multipoles with the moment symmetry axis tilted 
into the equatorial plane, i.e. perpendicular to the stellar ro- 
tation axis. 

As a test of the above Figure |4] shows simulated mag- 
netic surface maps for the four stars shown in Figure [T] The 
polar strength of the dipole and octupole components, the 
phase they are tilted towards, and their tilts relative to the 
stellar rotation axis match those listed in Table [T] For com- 
parison we also show in Figure|4]the observationally derived 
magnetic surface maps in the same Cartesian (i.e. latitude vs 
longitude) grid format. The £ number of a multipole deter- 
mines the number of polarity changes in the surface field if 
one moves along a line of constant longitude from the north 
to the south pole of the star (or vice versa). The domination 
of the dipole component (i = 1) in the A A Tau map and the 
octupole component (i=3) in the TW Hya map is clear. 

In order to produce field extrapolations some assump- 
tion has to be made about the form of the magnetic field 
in the ^hemisphere of the star that always remains hidden 
from view to an observer. In Figured] we have assumed that 
the field is dominantly antisymmetric with respect to the 
stellar midplane, and thus the odd i number modes domi- 
nate. If we assume that the even i number modes dominate 
then the resulting large scale field topologies would result in 
magnetospheric accretion occurring primarily on to lower, 
equatorial, latitudes. As this is typically not observed, and in 
order to ensure that material from the disk is able to accrete 
into the high latitude hot spots that are observed, we favor 
the antisymmetric field modes and assign the field in the 
hidden hemisphere accordingly. We note that this assump- 
tion has been examined in detail by Johnstone et al (2010) 
who conclude that for accreting T Tauri stars, the field in 
the visible hemisphere that is reconstructed in the magnetic 
maps remains largely unaffected by the assumptions of the 
form of the magnetic field in the hidden hemisphere. 

The effect of higher order field modes is apparent in 
the observationally derived surface magnetic maps shown 
in Figure m particularly for V2129 Oph and BP Tau. The 
high latitude single polarity magnetic spots are stronger and 
of different shape in the observationally derived maps com- 
pared to the simulated maps constructed from dipole and oc- 
tupole field components only. This is caused by the presence 




Fig. 5 Cartoon illustrating the definition of the height h 
and the width w of a. closed field line loop with foot points 
at colatitudes and longitudes (6>i, (/)i) and (6>2, (/>2). 

of higher order i numbers in the observed maps. Likewise 
there is more azimuthal structure apparent in the observed 
maps; in particular within the mid-latitude negative field 
band on both V2129 Oph and BP Tau. This is caused by 
the presence of non-axisymmetric field components (those 
with m ^ 0) in the observed maps (further details can be 
found in the papers where the magnetic maps have been 
pubHshed; Donati et al 2011a for V2129 Oph; Donati et al 
2008 for BP Tau; Donati et al 2010b for AA Tau; and Donati 
et al 201 lb for TW Hya). The maps derived from Zeeman- 
Doppler imaging thus contain more information than can 
be captured by the simple model consisting of tilted dipole 
plus tilted octupole field components only. None-the-less 
it is apparent from Figure (H that the simple dipole-octpole 
model provides excellent overall agreement on their large 
scale field topologies; it is such dipole-octupole models that 
have been used as the input conditions in the latest gener- 
ation of 3D MHD simulations of the star-disk interaction 
(Romano va et al 2011; Long et al 2011). In the following 
section we compare the large scale field topology derived 
via field extrapolation from the two sets of surface maps 
presented in Figure (H 



3.2.2 Comparison with field extrapolations 

We have carried out field extrapolations from the magnetic 
surface maps derived from Zeeman-Doppler imaging and 
those constructed assuming only tilted dipole and tilted oc- 
tupole field components shown in Figured This allows for 
a quantitative comparison between the 3D magnetospheric 
topologies derived from both sets of maps. To compare the 
field topologies we calculate the height and width of the 
closed field line loops, as defined in Figure \5\ The height 
h is the maximum height the loop reaches above the stellar 
surface, and the width w is the distance between the loop 
foot points as measured along the segment of the great cir- 
cle that passes through the foot points on the stellar surface. 
Details of how these quantities are calculated can be found 
in Gregory et al (2008). 
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Fig. 4 (left) Observationally derived magnetic maps assuming antisymmetric field configurations with respect to the 
stellar midplane for (top to bottom) AA Tau in January 2009, BP Tau in February 2006, V2129 Oph in July 2009 and 
TW Hya in March 2010. (right) simulated magnetic maps constructed using tilted dipole and tilted octupole components 
with the same polar field strength, tilts, and phases of tilt as determined from the observationally derived magnetic maps as 
listed in Table [T] 
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Figure [6] shows the comparison between the field ex- 
trapolations from the maps in Figure |4l For simplicity we 
have set the source surface to be at approximately the equa- 
torial corotation radius for each star. We remind readers 
that the source surface denotes the maximum extent of the 
closed field, and thus a field line cannot have a height greater 
than the radius of the source surface; in this case corotation. 
Larger scale field lines are likely torn open due to the in- 
teraction with the disk (e.g. Matt & Pudritz 2005), as has 
also been found in MHD simulations (e.g. Romano va et al 
2011). The maximum width that a field line could have is 
less than tt (loops with a larger width would be wider 
than the star, which is clearly unphysical). On each plot in 
Figures the solid red line shows the relationship between h 
and w for a pure dipole field, and the solid green lines for 
a pure octupole field. There are two such lines for the oc- 
tupole case as for an axisymmetric octupole there are three 
rings of closed field about the star, and the higher latitude 
rings are of different width compared to the equatorial ring 
(this is due to the mathematical properties of the associated 
Legendre functions, see Gregory 2011). 

It is clear from Figure [6l that there is good agreement 
between the 3D field structures obtained via field extrapola- 
tion from the observationally derived (left column), and the 
simulated dipole-octupole (right column), magnetic maps. 
Better agreement is however obtained for A A Tau (top row) 
and TW Hya (bottom row). These are the stars which have 
one of the dipole or octupole components clearly dominant 
over the other; AA Tau hosting a dominantly dipolar mag- 
netic field, and TW Hya a dominantly octupolar one (see Ta- 
ble [T]). The agreement is poorer, although not significantly 
so, for BP Tau and V2129 Oph (second and third rows of 
Figure O, stars which have both strong dipole and octupole 
field components; although for V2129 Oph/BP Tau the oc- 
tupole/dipole component contains the bulk of the magnetic 
energy (Donati et al 2011a, 2008). For V2129 Oph in par- 
ticular there are many field lines with large widths 2 R^) 
but small heights (~ I R^) that are apparent in the field ex- 
trapolation from the observationally derived map (see Fig- 
ure|4]first column third row) but which are missing from the 
extrapolation from the simulated map containing the tilted 
dipole and tilted octupole field components only (see Fig- 
ure |4] second column third row). These field lines connect 
the strong negative field regions evident in the mid-latitude 
northern hemisphere to equivalent positive field regions in 
the other hemisphere. This azimuthal field structure which 
arises due to the presence of non-axisymmetric (m ^ 0) 
field modes can be seen in the observationally derived map. 
As the simulated maps neglect the m / field modes for 
simplicity, a group of field lines is missed in the field extrap- 
olation from the V2129 Oph simulated map, which explains 
the difference seen in Figure [6l Given that these field lines 
are limited in height their presence does not significantly af- 
fect the structure of the larger scale field which is interacting 
with the disk, and is thus unlikely to have any effect on the 



results of MHD modeling which use dipole-octupole fields 
as inputs (see section (4]). 

AA Tau and TW Hya have field configurations where 
the dipole and octupole components are close to an anti- 
parallel configuration, whereas for BP Tau and V2129 Oph 
they are close to a parallel configuration (at least at the 
epochs listed in Table [T]). The differing field configurations 
are clearly apparent in the global field structure and in the 
plots of field line height versus width, see Figure [6j The 
largest field lines, those with the greatest height, for AA 
Tau and TW Hya have smaller width at the same height 
compared to a pure dipole, i.e. the black points in Figure 
[6] for these stars lie to the left of the red line that represents 
the dipole. The opposite is the case for BP Tau and V2129 
Oph where the field lines with the greatest height are wider 
than those of a dipole (as already discussed in Gregory et al 
2008) i.e. the black points in Figure [6l lie to the right of the 
red line. 

In all cases the magnitude of this effect, and the depar- 
ture from a dipole field, is greater for larger values of the 
ratio of the polar field strength of the octupole to the dipole 
component Boat/ B dip. This can be understood as follows. 
The dipole and octupole components for each of the four 
stars listed in Table [T] are tilted by different amounts to- 
wards different rotation phases. However, as the tilts are 
small, let's assume that V2129 Oph and BP Tau have dipole- 
octupole fields where the moments are both aligned with 
the stellar rotation axis and parallel to one another, /3dip = 
f^oct = 0° (main positive poles coinciding with the visi- 
ble rotation pole). Likewise, let's assume that AA Tau and 
TW Hya have dipole-octupole fields where the moments 
are both aligned with the stellar rotation axis but are anti- 
parallel to one another, Pdip = 0°, Poet = 180° (pos- 
itive/main negative pole of the dipole/octupole coinciding 
with the visible rotation pole; TW Hya is close to this con- 
figuration) or Pdip = 180°, Poet = 0° (negative/main pos- 
itive pole of the dipole/octupole coinciding with the visi- 
ble rotation pole; AA Tau is close to this configuration). 
The field components for the simplified cases of aligned 
and anti-aligned dipole-octupole fields can be derived from 
equations dMl]) using ([6l|8]), 

Br = Bdip (^"""^ ^ f^dip + 

1 / R \^ 

-Boot {~^] c^s^ ^ cos^ Poet - 3) COS 6 cos Poet (12) 

1 R \^ 

Be = -Bdip [~^] ^^^^ Pdip + 

3 / R \^ 

-Boot (~^] ^O^^ ^ ^^^^ f^oet - 1) sin 6> cos Poet (13) 

where 5^ = as the fields under consideration are ax- 
isymmetric, and the first/second terms in each component 
are the contributions from the dipole/octupole part of the 
field. Using these field components we can solve the differ- 
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Fig. 6 Height h versus width w plots for the field lines obtained via field extrapolation from the magnetic maps in 
Figure m (black points). Plots corresponding to the observationally derived maps are shown on the left, and those from the 
simulated maps assuming only tilted dipole plus tilted octupole field components are shown on the right for, from top to 
bottom, AA Tau, BP Tau, V2129 Oph and TW Hya. On each plot the red line shows the height- width relation for an axial 
dipole field (£ = 1, m = 0), and the green lines the equivalent for an axial octupole field (^ = 3, m = 0). 
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ential equation describing the path of the field lines, 

^ = ^, (14) 

dr rdO ' 

and determine the structure of the magnetic field external 
to the star. By path of the field lines we mean a function of 
the form r = r{6) that describes the field lines in spherical 
coordinates. For example, for an aligned dipole (set Boot = 
and Pdip = 0° in equations ([T2l) and ([T3]) ) the result is 
r = const, sin^ where the constant is the maximum radial 
extent of the dipole loops in the stellar midplane (a result 
commonly found in the literature e.g. Gregory 201 1). 

For the dipole-octupole fields considered here, and in 
fact for any arbitrary combination of axial £ number multi- 
poles, equation (O has an analytic solution; however, the 
mathematics is just too miserable to detail here. Equation 
([T4l) , with equations ([T2l) and ([T3]) , can also be solved nu- 
merically. Figure|7]shows the axisymmetric dipole-octupole 
magnetic fields for two cases (i) for dipole and octupole 
components that are aligned with the rotation axis and par- 
allel to one another with (3 dip = Poet = 0° (a config- 
uration similar to that of V2129 Oph and BP Tau), with 
Boot/ B dip = 4.0 to more clearly emphasize the difference 
between this field and a dipole; and (ii) for dipole and oc- 
tupole components that are aligned with the rotation axis 
and anti-parallel with (3 dip = 0° and Poet = 180° (a config- 
uration similar to that of TW Hya) again with Boat / Bdip = 
4.0. On each plot the red lines show some dipole field lines 
passing through the stellar midplane at the same position 
as some of the dipole-octupole field lines (black lines). The 
field topology depends only on the ratio Boot / Bdip rather 
than on the polar strengths of the individual field compo- 
nents. 

The field line plots in Figure [7] provide an immediate 
explanation for the differences found between the extrap- 
olated fields (the black points) and the dipole fields (red 
lines) in Figure [6l For stars like V2129 Oph and BP Tau 
where the dipole and octupole moments are roughly paral- 
lel, the large scale field lines are squeezed by the octupole 
field line loops as we move along the large scale field lines 
from the stellar midplane towards the stellar surface. This 
can seen in Figure [7] (top panel) by comparing the dipole- 
octupole (black) and the dipole (red) field lines. In such 
cases the large scale field lines of the dipole-octupole fields 
have larger width at the stellar surface than dipole field lines 
of the same height due to the squeezing by the more com- 
plex surface field regions. This effect becomes greater the 
larger the ratio of Boat/ B dip, i.e. the greater the influence 
of the octupole over the dipole component. This be seen in 
Figure [5] where, for a field line of a fixed height, we com- 
pare the width of the dipole-octupole field line relative to 
that of a dipole field line as we vary the ratio Boot/ B dip- 
As V2129 Oph has a larger Boot/ B dip ratio than BP Tau 
this squeezing effect is greater for V2129 Oph, and thus the 
large scale field lines of a given height are even wider than 
those of a dipole for V2129 Oph than they are for BP Tau 
(see Figures second and third rows and Figured]). For field 




x/R* 

Fig. 7 Field lines (black) for an axisymmetric dipole- 
octupole field where the dipole and octupole components 
are parallel with (3 dip = Poet = 0° (top panel), and where 
the dipole and octupole components are anti-parallel with 
Pdip = 0° and Poet = 180° (bottom panel). Boct/Bdip = 
4.0 in both cases. For comparison some dipole field lines 
are shown in red. The larger scale dipole-octupole field lines 
are squeezed by the more complex surface field. As the field 
lines are rotationally symmetric about the z-axis and reflec- 
tionally symmetric in the x-axis only one quadrant is shown. 

topologies like AA Tau and TW Hya, where the dipole and 
octupole moments are close to anti-parallel, the large scale 
field lines are distorted in a different way close to the stel- 
lar surface, see Figure [7] (bottom panel). In such cases the 
large scale field lines have smaller width at the stellar sur- 
face than dipole field lines of the same height, explaining 
the difference found in Figure [6] for AA Tau and TW Hya 
(first and fourth rows) between the extrapolated and dipole 
fields. As TW Hya has the largest Boct/Bdip value the ef- 
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Fig. 8 The variation in the ratio of the width of a large 
closed field line loop for an axisymmetric dipole-octupole 
magnetic field Wdipoie-octupoie relative to the width of a 
dipole loop Wdipoie of the same height as a function of the 
ratio of the polar field strength of the octupole to the dipole 
Boot I Bdip- For comparative purposes a height of h = AR^ 
has been chosen, which corresponds to a field line threading 
a disk at a radius of 5 i?*, a typical disk truncation radius. 
The width and height are defined in Figure [H The solid red 
line shows the case where the dipole and octupole compo- 
nents are parallel; in this senario, relevant to BP Tau and 
V2129 Oph, the large scale field lines are wider than the 
pure dipole case. The dashed blue line shows the case where 
the dipole and octupole and anti-parallel; in this senario, rel- 
evant to AA Tau and TW Hya, the large scale field lines are 
of smaller width than the pure dipole case. 



feet is greatest for this star, but is (almost) negligible for 
the dominantly dipolar magnetic field of AA Tau (the field 
line plot using the AA Tau Boot / Bdip ratio shows very little 
difference from a pure dipole). The results of this section 
demonstrate how an analytic model can be used to further 
our understanding of the results obtained from numerical 
(field extrapolation) models. 

The effect of the large scale field structure being dis- 
torted close to the stellar surface by the complex surface 
field regions is important for models of accretion flow along 
the field lines on to the stellar surface (Adams & Gregory 
2011). Not only is the hot spot location a sensitive func- 
tion of the magnetic field topology (e.g. Gregory et al 2005, 
2006a; Mohanty & Shu 2008), but the pre-shock density of 
gas can be up to an order of magnitude larger than that pre- 
dicted from dipolar accretion models (Gregory et al 2007). 
This is straightforward to understand physically if we as- 
sume that accretion occurs from the stellar midplane (e.g. 
imagine a flat disk in the equatorial plane of Figure [7]). As- 
suming mass conservation then the mass flux through the 



cross sectional area A of the accretion column is constant, 
i.e. pvA = const, where p and v are the flow density and 
velocity, and the constant is the mass accretion rate. The in- 
fall speed for gas coupled to the magnetic field is roughly 
the same for both accretion along dipole, and non-dipolar 
field lines, since the infall velocity is essentially free-fall 
and determined by the depth of the gravitational potential 
well of the star. However, the cross-sectional area of the 
dipole-octupole field lines becomes smaller than that of the 
dipole case (see Figure [7]), and therefore to ensure that the 
mass flux remains constant, the flow density must increase 
by roughly the same factor that the cross sectional area re- 
duces. 

3.2.3 Longitudinal field component 

When high-resolution Stokes V (circular polarization) data 
is available detailed modeling of the profiles and their rota- 
tional modulation provides the best information about stel- 
lar magnetic field topologies, and allows magnetic maps to 
be constructed. However, estimates of the stellar disk av- 
eraged longitudinal field component (also called the trans- 
verse or line-of- sight field component) Bion and how this 
varies with rotation phase ^ is commonly found in the lit- 
erature. The longitudinal field component is an integrated 
quantity over the visible disk of the star, and thus is often 
a poor diagnostic of the surface field topology due to the 
flux cancellation effect. If the stellar surface was covered in 
equal amounts of positive and negative field the measured 
stellar disk averaged longitudinal field component would be 
zero, yet in such cases the Stokes V profile itself may show a 
clear Zeeman signature (e.g. see the longitudinal field curve 
from the LSD average photospheric absorption line and the 
Stokes V profile for V2129 Oph around phase 0.78 in July 
2009; Donati et al 2011a). 

Large scale magnetic topologies have been probed us- 
ing the oblique rotator model, where a dipole field is as- 
sumed to be tilted by an angle Pdip relative to the stellar 
rotation axis with a stellar inclination of i. The longitudi- 
nal field component thus varies with rotation phase, reach- 
ing a maximum when the star is viewed at a rotation phase 
^ that corresponds to the phase that the dipole component 
is tilted towards ^dip (assuming it is the positive pole of 
the dipole in the visible hemisphere). Here we extend our 
analytic work by deriving an expression for the longitudi- 
nal magnetic field component for dipole-octupole magnetic 
fields. Even for such simple field structures, and where the 
octupole component is stronger than the dipole component, 
we demonstrate below that the phase modulated longitudi- 
nal field curve is dominated by the dipole component (al- 
though the polar field strength and tilt of this component is 
also poorly constrained by the longitudinal field curve once 
cool spots, neglected in this simple model, are accounted 
for - see below), and therefore provides poor constraints on 
the stellar field topology. Detailed modeling of rotationally 
modulated Zeeman signatures, when high resolution spec- 
tropolarimetric data is available, provides better and more 
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detailed constraints on stellar magnetic topologies (see Do- 
nati & Landstreet 2009 for further details on the limitations 
of the longitudinal magnetic field component). 

Although the expression for the longitudinal field com- 
ponent for an inclined dipole magnetosphere is commonly 
found in the literature (e.g. Preston 1967), we cannot find 
explicit expressions for a tilted octupole magnetosphere, al- 
though we note that Bagnulo et al (1996) did use a combi- 
nation of a tilted dipole, a tilted quadrupole and a tilted oc- 
tupole in their work. Assuming a linear limb darkening law 
the stellar disk averaged longitudinal magnetic field com- 
ponent for a dipole-octupole field as a function of rotation 
phase is (see Appendix |B] for a derivation of this result), 

{cos j3dip COS i + sin Pdip sin i COS [27r($ - ^dip)]) 

|5(cos /3oct cosz + sin /3oct sin i cos [27r($ - ^oct)])^ - 

3(cos f3oct cos i + sin f3oct sin i cos [27r(<l> - ^oct)]) } , (15) 

where B^ip is the polar strength of the dipole component, 
Pdip the tilt of the positive pole of the dipole relative to 
the visible stellar rotation pole, ^dip is the rotation phase 
that the dipole moment is tilted towards, with Boot-, Poet 
and ^oct the same respective quantities for the octupole 
field component. $ is the rotation phase, and u is the lin- 
ear limb darkening coefficient. The stellar inclination is z, 
such that the observer's reference frame {xobs^Vobsj ^obs) 
has z-axis parallel to the line of sight i.e. we can transform 
from the stellar reference frame {x^y^z) to the observer's 
frame {xobs^Vobs^ Zobs)^ which share the stellar center as 
the coordinate origin, by rotating the stellar frame counter- 
clockwise by an angle i when looking down the ?/-axis to- 
wards the origin. 

In Figure [9] we illustrate the longitudinal field compo- 
nent as a function of rotation phase for a star with a dipole- 
octupole field with the same parameters as BP Tau listed 
in Table [T] It is clear that the longitudinal field curve is 
dominated by the dipole component of the magnetic field 
even though the octupole field component is stronger. Ob- 
servationally the octupole component may easily be missed 
by analyzing the longitudinal field curve alone once a lim- 
ited phase coverage and errors are accounted for, despite 
being easily recoverable if a full analysis is made of the 
rotational modulation of the Stoke V signal, as is done in 
Zeeman-Doppler imaging and in the construction of stel- 
lar magnetic maps. The situation is worse for higher order 
field components. Furthermore, even though the longitudi- 
nal field curve is only weakly sensitive to the high order 
field components in the simple model presented here, it can- 
not be used as a probe of the large-scale dipole component. 
In this simple analytic dipole-octupole field model we have 
not accounted for the presence of surface cool (dark) spots. 



0.4 I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




I I I I I I I I I I I I I I I I I 

0.2 0.4 0.6 0.8 1 
rotation phase (4>) 

Fig. 9 The stellar disk averaged longitudinal magnetic 
field component Bion as a function of rotation phase ^ for 
a star with a tilted dipole plus tilted octupole magnetic field 
with tilts, phases of tilt, and polar field strength of the com- 
ponents matching those determined from the observation- 
ally derived magnetic map of BP Tau, see Table [TJ Bion 
(black line) is calculated from equation (fTSl) assuming a lin- 
ear limb darkening coefficient of u = 0.6 and a stellar in- 
clination of 45°. The red/blue line shows the contribution to 
Bion from the dipole/octupole component of the field. This 
simple model ignores surface dark spots, which should be 
accounted for, see the text for details. 

Brightness maps are produced as part of the magnetic imag- 
ing process, and for many accreting T Tauri stars the main 
magnetic pole often coincides with large dark spots (e.g. for 
BP Tau; Donati et al 2008b). Dark spots coinciding with the 
main magnetic pole will significantly change the respective 
contributions of the dipole and octupole field components 
to the longitudinal field curve, with the octupole producing 
a comparably larger contribution. If the cool spot is off- set 
from the main magnetic pole then the tilt of the dipole com- 
ponent /3dip that can be derived by fitting equation (fTSl) to an 
observed longitudinal field curve will likely be significantly 
off. 

If the longitudinal field curve constructed from a LSD 
average photospheric line was probing the large scale dipole 
component of a stellar magnetic field, then the variations 
of Bion with phase would be mostly sinusoidal; observa- 
tionally this is not usually the case. For V2129 Oph the de- 
parture from a simple sinusoidal longitudinal field curve is 
obvious, see the plot in figure 2 (bottom left panel) in Do- 
nati et al (201 la). For AA Tau the longitudinal field compo- 
nent measured from the LSD average photospheric absorp- 
tion line is almost featureless, and yet detailed modeling of 
the full Stokes V profile reveals a strong dipole component 
of polar strength ~ 2 kG. Longitudinal field components 
therefore provide poor observational constrains on T Tauri 
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magnetic field topologies even for the dipole component it- 
self. 

Although strong and rotationally modulated longitudi- 
nal field components are often measured in accretion related 
emission lines (e.g. Valenti & Johns-Krull 2004; Donati et al 
2008b), there is another reason why fitting an oblique dipole 
rotator model to the longitudinal field curve may yield po- 
lar field strengths for the dipole component that are signifi- 
cantly discrepant from the true value. If we use V2129 Oph 
in July 2009 as a specific example, the dipole component 
recovered from a spherical harmonic decomposition of the 
magnetic map suggests a dipole component of ~ 1 kG, but 
the bulk of the material arrives at the star into a single pos- 
itive field spot of ^ 4kG (Donati et al 2011a). Thus the 
analysis of the longitudinal field component measured in 
accretion related emission lines which are forming close to 
the accretion shock would yield discrepant values for the 
dipole component (on V2129 Oph higher order field com- 
ponents contribute to the strong field region at the base of 
the accretion funnel). Indeed Valenti & Johns-Krull (2004) 
find that their longitudinal field components measured in 
the Hel 5 876 A emission line are rotationally modulated, 
but well fitted by a simple model where the accreting ma- 
terial lands in a single polarity spot on the stellar surface. 
This is found in Zeeman-Doppler imaging studies for many 
accreting T Tauri stars (e.g. Donati et al 2008b; Donati et 
al 2011a). Analysis of the longitudinal field curve alone 
(whether derived from accretion related emission lines or 
from a LSD average photospheric absorption line) therefore 
provides poor constrains on the field topology and the po- 
lar strength of the dipole component for accreting T Tauri 
stars. A full analysis of variations in the Stokes V profile for 
both accretion related emission lines, and photospheric ab- 
sorption lines, is clearly desirable whenever high resolution 
spectra are available to provide a consistent description of 
the large-scale magnetic field topology. 

4 Discussion and conclusions 

In this paper we have discussed the numerical field extrapo- 
lation and analytic models employed over the past few years 
to construct 3D models of the magnetospheres of accret- 
ing pre-main sequence stars. Each of the methods, analytic, 
field extrapolation, and also magnetohydrodynamic which 
we discuss below, have their own advantages, disadvantages 
and limitations. By writing out the field components for 
magnetospheres that consist of several high order field com- 
ponents analytically, as in section 13.21 results that may be 
missed by numerical models are revealed. As one example, 
if we consider the simplified case of a dipole plus octupole 
magnetic field, where the dipole and octupole moments are 
both aligned, then it becomes clear from equations ([T2l) and 
([T3]) that a magnetic null point exists where Br = Bq = 
at a radius of 



^nuii I R* is the radius of circle in a plane tilted by /3 (the 
tilts of the dipole and octupole moments) towards the phase 
that the multipole moments are tilted towards. Essentially 
this radius defines the point at where the dipole compo- 
nents begins to dominate the octupole component as we 
move away from the stellar surface towards larger radii. 
As the polar strength of the octupole relative to the polar 
strength of the dipole. Boat/ B dip, is increased the extent 
of the region where the octupole component dominates in- 
creases and "smaller scale" field lines extending to greater 
height are found in the field extrapolations; this can be seen 
from the plots in Figure [6] for stars with larger Boat / Bdip 
ratios (see Table [T]). 

The major disadvantage of the analytic approach is the 
inclusion of a large amount of £ number multipoles rapidly 
becomes cumbersome and such models are typically (but 
not exclusively) limited to axisymmetric fields (e.g. Mo- 
hanty & Shu 2008; Gregory et al 2010; Gregory 201 1). Fur- 
thermore, when the multipole moments are tilted by differ- 
ent amounts relative to the stellar rotation axis, and/or to- 
wards different rotation phases, the field structure (i.e. the 
paths of the field lines exterior to the star) must be derived 
numerically from 5^ /dr = Bo/{rdO) = ^^/(r sin6>d0). 
Inclusion of non-axisymmetric field modes would further 
add to the complexity of the analytic work. In such cases 
models of the magnetospheres of stars with complex mag- 
netic fields are best handled numerically. 

Numerical field extrapolations from magnetic surface 
maps have the advantage over the analytic approach that as 
many £ and m modes that can be resolved observationally 
can be included. Field extrapolation is currently the only 
method where magnetic fields with an observed degree of 
complexity can be handled and incorporated into models 
of the large-scale magnetosphere. Potential field models are 
computationally simple to carry out and only require desk- 
top computing resources. They also produce unique field 
topologies (Aly 1987). However the field topologies ob- 
tained via potential field extrapolation are static. Time de- 
pendent effects, such as the distortion of the large scale 
magnetic field due to the interaction with the circumstellar 
disk, cannot be modeled. Other non-potential effects, such 
as flux emergence or surface differential rotation that has 
been measured for some accreting T Tauri stars (e.g. Donati 
et al 2010a, 201 la), also cannot be included. 

Three-dimensional MHD models of the star-disk inter- 
action with non-dipolar magnetic fields have also been de- 
veloped by Long et al (2007, 2008, 2009, 2011) and Ro- 
mano va et al (2011). The main advantage of the MHD ap- 
proach is the ability to handle the time evolution of both 
the large-scale field, and the magnetospheric mass accretion 
process, that the potential field extrapolation models cannot. 
Unfortunately 3D MHD simulations require supercomputer 
resources; such models are therefore computationally, and 
fiscally, expensive. Likewise, magnetic maps have not yet 
been incorporated into MHD simulations as grid resolution 
issues arise. A small grid resolution close to the star is re- 
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quired in order to include high order and non-axisymmetric 
field components, but at the same time, a large grid able to 
reach out as far as the inner disk is required. Instead the lat- 
est generation of 3D MHD models, presented by Romanova 
et al (2011) and Long et al (2011), begin by assuming that 
the star has a potential magnetic field consisting of a tilted 
dipole plus a tilted octupole component with field strengths 
and tilts being adopted from the published magnetic map 
of the star of interest. Such dipole-octupole field structures 
are then allowed to evolve due to the interaction with the 
disk. The largest scale field lines are torn open and quickly 
twisted around the rotation axis of the star, see e.g. figure 6 
of Romanova et al (201 1), however, the field interior to this 
that carries gas in columns on to the star retains its initial 
potential structure. 

MHD models that incorporate the magnetic maps de- 
rived from Zeeman-Doppler imaging remain to be carried 
out in the future. Likewise a major new avenue for research 
is the inclusion of stellar surface differential rotation, flux 
emergence and meridional flows into MHD models. It may 
turn out that the retention of the potential field structure 
found in current simulations is no longer valid once stellar 
surface transport effects have been accounted for. Likewise 
it is unclear whether T Tauri stars undergo magnetic cycles, 
although there is tentative evidence that they do and that 
their large scale magnetospheres evolve over long (possibly 
several years) timescales (e.g. Donati et al 2011a). 

The veracity of the 3D field structures obtained via field 
extrapolation from the magnetic maps must continue to be 
examined in future. On going large multi- wavelength pro- 
grams are useful here (e.g. Argiroffi et al 201 1; Donati et al 
201 Ic). Based on a field extrapolation model of V2129 Oph 
Jardine et al (2008), using a simple accretion flow model 
presented in Gregory et al (2007), predicted that a clear 
signature of accretion related X-ray emission would be de- 
tected, which was subsequently found (Argiroffi et al 201 1). 
Never-the-less, further tests are underway to determine if 
the derived field topologies genuinely capture the true mag- 
netospheric structure of T Tauri stars. 
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A Tilted multipole moments in 3D 



stellar frame (x,y,z) Q- Oi multipole frame (x',y',z') \k = nz' 
z 




Fig. Al A given multipole moment symmetry axis /j. is 
assumed to be tilted by an angle (3, relative to the stellar 
rotation axis, towards a longitude ip, where is related to 
the rotation phase via = (1 — phase) x 360°, assuming 
the xz-plane corresponds to a rotation phase of zero. The 
stellar reference frame (x^y^z) is defined such that ft = 
ftz, and the multipole frame (x', y' ^ z') such that fi = /iz'. 
We can transform from the stellar coordinate frame to the 
multipole frame by 1 . rotating the stellar coordinate system 
by i/j around the z-axis, followed by 2. a rotation about the 
new y-Sixis by an angle (3 as illustrated. 



with the 2;'-axis (/j, = /xz^. For an arbitrary ^-number multipole 
Gregory et al. (2010) showed that the magnitude of the multipole 
moment is related to the polar strength of the multipole via, 



(Al) 



We assume that the xz-plane corresponds to a rotation phase of 
zero. The multipole moment is assumed to be tilted by an angle /3, 
such that tt • [l — cos /3, towards a phase determined by another 
angle V^, the longitude in the stellar reference frame, and where 
Q. and ji are unit vectors along Q, and ^ respectively. In other 
words 13 — 45° and ^ = 90°, represents a multipole moment 
tilted by 45° from the stellar rotation axis, towards rotation phase 
0.75 i.e. = (1 — phase) x 360°. Thus in order to transform 
from the coordinate system that denotes the frame of the star to 
that which denotes the frame of the multipole we first rotate coun- 
terclockwise about the z-axis (when looking towards the origin) 
by an angle ^, and then counterclockwise by an angle 13 (when 
looking towards the origin) about the new ^/-axis, as illustrated in 
Figure IaT] (note that we are rotating the entire coordinate system, 
and not simply a vector within the same fixed coordinate system). 
The Cartesian components of the multipole moment in the stellar 
reference frame, where ^ — /j.x^ + /J^yY -\- /J^zZ, are 




cos p sin /3 

1 
— sin P cos P 




(A2) 



We consider two Cartesian coordinate systems (i) the stellar ref- 
erence frame (x^y^z) where the stellar rotation axis is aligned 
with the z-axis (ft — Qz) and (ii) the multipole reference frame 
(x^ y' , z) where the multipole moment symmetry axis is aligned 



The vector components of /j, can then be converted from the Carte- 
sian to a spherical coordinate system (r, 0, (j)) where r = corre- 
sponds to the center of the star, = corresponds to the stellar 
xz-plane and is the co-latitude in the stellar reference frame. 
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Noting that /it = /Xa^x + /i^y + /i^z = /x^r -\- /j^qO -\- /j.cf)(t>, then 
sin ^ sin cos^ \ / \ 
cos ^ sin — sin ^ /i^ (A3) 
cos (/) J \ fiz J 

The spherical components of the tilted multipole moment in the 
stellar reference frame are thus given by equations (I6ll8] |. 

B Stellar disk averaged longitudinal 
magnetic field component 

Goossens (1979) derived an expression for the stellar disk aver- 
aged longitudinal magnetic field component Bion arising from a 
multipolar stellar magnetic field, see his equation (20) with the ' W 
terms for a linear limb darkening law from his equations (25a-d) - 
see also Stibbs (1950) and Schwarzschild (1950). It is applicable 
to the case where all of the multipole moments are aligned with 
each other, and which are viewed at an angle a to the magnetic 
axis. The first term in Goossens' equation (20) is the dipole term, 
the second term is the quadrupole term, the third term represents 
the contribution from the higher order even £ number multipoles 
(i = 4 [hexadecapole], 6, 8 . . .), and the fourth term is the con- 
tribution from the higher order odd £ number multipoles {£ — 3 
[octupole], 5 [dotriacontapole], 7 . . .). Goossens then goes on to 
derive the stellar disk averaged longitudinal field component for a 
dipole-quadrupole magnetic field, but this can be extended to any 
arbitrary combination of high order £ number multipoles. 

By comparing the magnetostatic potential in Goossens (1979) 
[his equation (4)] with that used in Gregory et al. (2010) [their 
equations (3.18) and (3.27)], the 'A' terms in Goossens (1979) 
equation (20) are related to the polar field strengths of the par- 
ticular multipoles via. 



An = Mf 



{£ + 1 



73^, pole 



where is the polar field strength of the ^th multipole, and 

we rewrite 5^'^°^^ = Bd^p, ^2'^°'^ = B^uad, and 5^'^°^^ = 
Boot for brevity. Thus, the 'A terms in Goossens equation (20) 
can be re- written as 

Al — -RlBdip 



A 



A2 
A2k 
2fc + l 



2 -R* Bquad 



1 j^2/c + 2^2/c,pole 



(2/c+l) 
1 



^2/c+3^2/c + l,pole 



{2k + 2) 

where the first term is the dipole term, the second the quadrupole 
term, the third term the higher order even £ number multipoles, and 
the fourth term the higher order odd £ number multipoles. Equation 
(20) of Goossens can then be re-expressed by replacing the 'A and 
'W terms, 

1 . . 



20 (3 - u) 
1 _ u 

3u 



~\~ . Bquad \ 

{3-u) 



P2(cos a) 



+ 



B 



2fc,pole 



■P2fc(cos a) 



^2/c + l,pole 



3(1 -u) 



P2fc+i(cosa), (Bl) 



where Pj (cos a) are the Legendre polynomials, a is the angle be- 
tween the magnetic axis (here it is assumed all the multipole mo- 
ments are aligned with each other) and the observer's line-of- sight, 
and 

_ fe+i (2/c + l)(2/c-3)!! 
- ^"^^ ¥{kT2)\ 



Dk = (-1 
with (2A;-1)!!: 



2fc(A; + 2)! 
(2k-l)(2k-3){2k- 



5)...5 x3 X 1. 



B.l Dipole-octupole longitudinal magnetic field 
component 

The dipole longitudinal field component is the first term of equa- 
tion dBlb . and derivations can be found in Schwarzschild (1950) 
and Stibbs (1950). The octupole magnetic field term is the /c = 1 
term in the final summation term of equation dBlb . Thus the disk 
averaged longitudinal magnetic field component for a pure oc- 
tupole is. 



B?, 



1 R in- 
16 (3- 



1) 



(5 cos^ a — 3 cos a) 



where a is the angle between the octupole moment and the ob- 
server's line-of- sight, a is related to the stellar inclination i and 
the angle between the stellar rotation pole and the multipole mo- 
ment [3 via. 



cos a = cos i cos /S + sin i sin /S cos c 



(B2) 



where is the azimuthal angle (the longitude) of the magnetic 
axis at some time t i.e. = 27r(l — $) where $ is the rotation 
phase (see the discussion in Goossens (1979) around his equations 
(23) and (24)). Equation (IB21 ) can be derived by considering the 
circle swept out by the magnetic axis rotating about the stellar ro- 
tation axis (it is just the spherical cosine law). If at time t = 0, 
i.e. rotation phase $ = 0, the magnetic axis is tilted towards some 
rotation phase $o then cos = cos (27r[l — + 27r[l — $o]) = 
cos (27r[$ - $o]), thus, 

cos a = cos i cos /3 + sin i sin /3 cos (27r — $o] ) • (B3) 

The stellar disk averaged longitudinal magnetic field component 
for an octupole is then, 

BLtW = ( J^) X 

|5(cos Poet cosi + sin Poet sin i cos [27r($ - ^oet)])^ 
—3{cos Poet cosi + sin Poet sin i COS [27r($ — $oct)])|. (B4) 

Goossens (1979) consider the dipole-quadrupole case where 
the dipole and quadrupole moments are tilted by different amounts 
towards different rotation phases (there is a sini term missing 
in the second term of his equation (35)). The Goossens work is 
straightforward to extend to dipole-octupole magnetic fields, with 
the dipole and octupule tilted towards different rotation phases and 
by different amounts relative to the stellar rotation axis. Taking the 
dipole and octupole components of equation dBll ) 

BlonW = — _ COS adip 

+ -^Boetj^ — ^ (5 cos^ ttoct - 3 COS aoct)(B5) 
16 (3 — u) 
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where adip and aoct are the angles between the dipole moment 
and the stellar rotation axis, and the octupole moment and the stel- 
lar rotation axis, respectively. In the notation adopted by Goossens, 
see their equation (35), A is the angle between the plane that con- 
tains the dipole moment and the stellar rotation axis, and the plane 
that contains the quadrupole moment (octupole moment in our 
case) and the stellar rotation axis. Thus, in the Goossens notation, 
at time t = if the dipole moment is titled towards longitude 00 
then the octupole (or quadrupole in the Goossens paper) moment 
is tilted towards longitude 00 + A; in other words we can write, 

cos adip — 

cos i cos Pdip + sin i sin /Sdip cos (27r[$ - ^dip]) (B6) 
cos aoct — 

cos i cos /3oct + sin i sin Poet cos (27r[$ - $oct]). (B7) 

Putting the equations in this section together gives equation ( fTSt , 
the general result for the dipole-octupole stellar disk averaged lon- 
gitudinal magnetic field component as a function of rotation phase. 
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